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Abstract 

We analytically calculate the contribution of transverse momentum con- 
servation to the azimuthal correlations that have been proposed as signals for 
possible local strong parity violation and recently measured in heavy ion col- 
lisions. These corrections are of the order of the inverse of the total final state 
particle multiplicity and thus are of the same order as the observed signal. 
The corrections contribute with the same sign to both like-sign and opposite- 
sign pair correlations. Their dependence on the momentum is in qualitative 
agreement with the measurements by the STAR collaboration, while the pseu- 
dorapidity dependence differs from the data. 
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1 Introduction 



Topological configurations occur generically in non-Abelian gauge theories and are 
known to be essential for understanding the vacuum structure and hadron proper- 
ties in Quantum Chromodynamics (QCD) (for reviews see e.g. p]). They have also 
been shown to play important roles in hot QCD matter (i.e. quark-gluon plasma) 
existing in the early universe and now created in heavy ion collisions [2J. Despite 
many indirect evidences, a direct experimental manifestation of the topological ef- 
fects has not been achieved and is therefore of great interest. One salient feature 
of the topological configurations is the V- and CP-odd effects they may induce. 
Based on that, it has been suggested [3j HI [6] to look for possible occurrence of 
V- and CP-odd domains with local strong parity violation for a direct detection 
of topological effects. Such domains may naturally arise in a heavy ion collision 
due to the so-called sphaleron transitions in the created hot QCD matter. In par- 
ticular, the so called Chiral Magnetic Effect (CME) predicts that in the presence 
of the strong external (electrodynamic) magnetic field at the early stage after a 
(non-central) collision, sphaleron transitions induce a separation of negatively and 
positively charged particles along the direction of the magnetic field which is perpen- 
dicular to the reaction plane defined by the impact parameter and the beam axis. 
Such an out-of-plane charge separation, however, varies its orientation from event to 
event, either parallel or anti-parallel to the magnetic field (depending whether the 
CME is caused by sphaleron or antisphaleron transition). As a result the expecta- 
tion value of any P-odd observable vanishes and only the variance of such observable 
may be detected, making the measurement of CME rather challenging. Recently 
the STAR collaboration has published in [7j measurements of charged particle az- 
imuthal correlations proposed in [8] as a CME signal, and found interesting patterns 
partly consistent with CME expectations. The STAR data have generated consid- 
erable interests and many subsequent works have appeared, proposing alternative 
explanations [91 [10J [EH 12], suggesting data interpretations and new observations 
[E2 QU HSl M\, and studying further consequences of CME [TH QH [T9] . 

We start with a discussion of the proposed CME signal as measured by STAR. 
In Ref. [8] it has been suggested that the CME may be indirectly approached by 
the measurement of the following two-particle correlation 

7 = (cos(0i + <p 2 - 2^rp)} 

= (cos((j) 1 -y R p)cos((P 2 -y R p))-{sm((j) 1 -y R p)sm((j) 2 -y R p)), (1) 

where ^ R p, and 02 denote the azimuthal angles of the reaction plane and pro- 
duced charged particles, respectively. Since this observable measures the difference 
between the in-plane and out-of-plane projected azimuthal correlations, it has been 
argued that this observable is particularly suited for revealing the CME signal which 
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is an out-of-plane charge separation. Specifically, the CME predicts 7 > for op- 
posite sign-pairs and 7 < for same-sign pairs. However, as has been pointed out 
in 0, [TOj [TTJ [TH [151 [16], in non-central collisions the presence of elliptic flow [20] 
already differentiates between in- and out-of-plane. As a consequence, essentially 
any two-particle correlation will contribute to the above observable even though the 
correlations' dynamical mechanism may generally be reaction-plane independent. 
These so-called "background" correlations, need to be well understood theoretically 
and possibly be determined by independent measurements^ The STAR publication 
[7] presents data for the correlator 7 together with the reaction plane independent 
correlator 



in the midrapidity region for both same- and opposite-sign pairs in AuAu and CuCu 
collisions at two energies ^/snn = 200 and 62 GeV. The data are encouraging: at 
first sight the results for 7 seem to be qualitatively consistent with the CME expec- 
tations. However, as shown in [U] when taking the correlator S into account as well, 
the interpretation of the data in term of the CME requires almost exact cancellation 
of the CME and all possible "background" correlations. Consequently in order to 
extract a possible signal for the CME, the understanding of these "background" 
correlations becomes crucial at this stage. 

One well-known possible source of azimuthal correlation is the conservation of 
transverse momentum, which has been qualitatively discussed in [11] and has been 
suggested to be a significant contribution to the measured observable 7. The argu- 
ment goes as follows. Consider for a moment all particles in the final state, charged 
and neutral over all phase-space. Next rewrite the correlator (TT]) as (for simplicity 
we set ^rp = here and in the rest of the paper) 



where i and j are summed over all particles. If we further assume that all particles 
have exactly the same magnitude of transverse momentum p t , the conservation of 
transverse momentum implies 



1 It is worth emphasizing that the contributions to 7 from elliptic flow induced correlations 
and from the CME have similar centrality trends: both the elliptic flow and the magnetic field, 
necessary for the CME, increase from central to peripheral collisions. 




(2) 




(3) 




(4) 
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and in consequence one obtains for sufficiently large N 



Here v-i is the elliptic flow coefficient measured for all produced particles and N is 
the total number of all produced particles (in full phase space). This contribution 
to the azimuthal correlations from transverse momentum conservation turns out to 
be at the order of the data measured by STAR, and therefore bears interest and 
importance. 

However, the above argument relies on two assumptions which are not realized in 
the actual measurement. First, STAR measures only the charged particles in a small 
pseudorapidity region \r]\ < 1, accounting only for about 15% of the total number 
of produced particles. Second, the magnitude of the transverse momentum is not 
a constant but rather distributed, more or less according to a thermal distribution. 
Therefore, a more realistic estimate for the contribution of transverse momentum 
conservation to the above correlation functions is required. The influence of momen- 
tum conservation on observables in heavy ion collisions has been discussed in the 
literature in the context of spectra [21] , elliptic flow [221 [21] and certain two-particle 
densities [23], and corrections of various importance have been established. In the 
present paper, we will address the effect of transverse momentum conservation on 
the correlation functions relevant for the potential measurement of the CME. To this 
end we derive the necessary formalism which allows us to quantify the azimuthal 
correlations due to transverse momentum conservation. 

Before going into details, let us discuss a few qualitative features which are to be 
expected from transverse momentum conservation, and which will be demonstrated 
in detail below. First, transverse momentum conservation introduces a back-to- 
back correlation for particle pairs, as they tend to balance each other in momentum. 
Second, the expected correction should scale inversely with total number of particles, 
as more particles provide more ways to balance the momentum and thus dilute the 
effect on two-particle correlations. Furthermore the correlation should be stronger 
in-plane than out-of-plane due to the presence of elliptic flow. As a result we expect 
that transverse momentum conservation results in a negative contribution to the 
observable 7, which increase with the strength of the elliptic flow, v-i- Finally, we 
note that transverse momentum conservation is "blind" to particle charge, leading 
to identical contributions to same-sign and opposite-sign pair-correlations. Because 
of these features, transverse momentum conservation alone cannot be expected as 
a full account for the observed charged particle azimuthal correlation patterns. It 
should be rather considered as an important background effect that contributes 
significantly and, therefore, necessitates quantitative studies for establishing any 
final interpretation of the data. 
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The paper is organized as follows. In the next Section, assuming only transverse 
momentum conservation, we present analytical calculations of the correlators 7 and 
5. In Section 3 we take into account the STAR acceptance and compare our results 
with the available data (integrated over transverse momentum). In Section 4 we 
discuss specific differential observables, such as transverse momentum and rapidity 
dependent correlations. Our conclusions are listed in the last section, where also 
some comments are included. 



2 General formulas 

Let us assume that there are N particles in total produced in a given heavy ion 
collision event, with individual momenta pi, ...,Pn- We denote the particles' trans- 
verse momenta by pi <t , ■ ■■,PN,t and their magnitudes by pi )t , ...,pN,t- The N particle 
density (normalized to unity) with enforced transverse momentum conservation 
can be written as 

_» s = 5 2 (p ht + ■■■ +PN,t) f(Pl)-f(PN) ^ 

f F S2 (Pht + ■■■ +PN,t) f(pi)-f(PN)d 3 Pi...d 3 p N ' 

where f(pi) is the normalized (J F f(p)d 3 p = 1) single particle distribution. Note that 
the above integrals are taken over the full phase space (denoted by U F" ) rather than 
the region where particles are actually measured. In Eq. (jBJ) we explicitly assume 
that all produced particles are governed by the same single particle distribution q 
We also ignore any other two-particle correlations, since in this paper we focus on 
the effects of transverse momentum conservation only. To calculate a two-particle 
correlator, such as (cos(0i + <p2))i we need the two-particle density which can be 
obtained from Eq. (jf)]) by integrating out N — 2 momenta 

h = f(Pi)f(p 2 ) )— f— • (7) 

f F * (El a) nil m)<m 

In order to perform the integrals in Eq. ((7j) we follow the techniques of Refs. 
2U and make use of the central limit theorem. The sum of M uncorrelated 



transverse momenta Yl^=iPi,t = K-t nas a Gaussian distribution if M is sufficiently 



2 This assumption is reasonably well satisfied in heavy ion collisions, where the final state par- 
ticles are mostly pions. 
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largcl that is 

G M (K t )= j 5 2 (j^ti^-^Uti [/(ft)^] 

J F 



2 



exp 



K 2 K 
ly x y 



Here x and y denote the two components of transverse momentum and 

(PI) F = l f(P)Pld 3 p, (PI) F = l f(p)p 2 yd 3 p, (9) 



where the integrations are over full phase space F. Using Eq. ([8]) we can express 
Eq. ((7j) in the following way 

G N - 2 {~Pl,t ~ P2,t) 



f2(pi,P2) = f(Pl)f(P2 



G N (0) 



- f(^)f(^)J^ exv f (P^+P^) 2 (Pi,y+P2,y) 2 ) m 

-nPi)HP2) N _ 2 ewy 2{n _ 2){p 2 )f 2(A r- 2 )<^>J- (10) 

Expandingin powers of 1/N and restricting ourselves to pairs with |pi )t + P2,t\ "C 



y/2N {p 2 ) U we obtain 

, 2 (p liX +p 2 ,x) 2 {Pl,y +P2,y) 2 \ m s 

Mn,P2) - /M/M (i + j, - 2N(A)f - 2N{p1)f J • (ID 

As already mentioned, the above correlation function does not distinguish between 
same- and opposite-sign pairs (or pairs involving neutral particles), as the transverse 
momentum conservation involves all particles without discrimination on the charge 
of particles. We also note that the two-particle density, f2(pi,P2) is maximum 
for back-to-back configurations, as can be easily seen from the dependence on the 
combinations —(j)i )X + p 2 , x ) 2 and ~(pi, y +P2, y ) 2 - 

Given the two-particle density Eq. ( ITlfTTT) we can proceed to evaluate various 
two-particle azimuthal correlations in any given kinematic region, for example the 
one introduced in Eq. (GQ) 

/ MPi,^) cos(0i + 2 )d 3 Pid 3 P2 

(COS(0 1 + 2 )> = ° ^W3^3- > ( 12 ) 

J j2{Pi,P2)d i pid i p2 
n 



3 We have checked by simple MC calculations that the central limit theorem applies for M > 10, 
see also [21) . This condition is very well satisfied in heavy ion collisions, even in the most peripheral 



ones. 
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In practice this is hardly any restriction, since for sufficiently large N only pairs with very 
large transverse momentum violate the condition \p\ t t + P2,t\ *C y/2N (pf) , which are strongly 
suppressed by the exponentially decreasing single particle distributions, f(pi)f(j>2)- 
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(cos(fa + fa)) = -—t J^T ' /= ,2 , (14) 



and analogously for (cos(0i — ^2))- Here we denote by Q the part of the phase space 
covered by the actual experiment. 

Finally we need the single particle distribution, for which we assume the following 
rather general form (^rp = 0): 

f(p^d 3 p=^^[l + 2v 2 (p t , V )cos(2fa}d 2 p t d V , (13) 

where v 2 is the pt and rj dependent elliptic flow coefficient (with r\ the pseudora- 
pidity). Taking Eq. ( ITT1) into account and performing elementary calculations we 
obtain our main result □: 

1 (p t ) 2 n 2v 2t n - h,F - ^2 ,f(^ 2 ,q) 
N (p 2 t ) F 1 - (V2 tF ) 

where we have introduced certain weighted moments of v 2 

_ 2 = (v2(j>t,v)Pt) = 1 9{Pt, rpvijpt, 7])p t d 2 p t dr) 
(Pt) J g(p t ,v)Ptd 2 Ptdr] 

and 

= 2 = (v2(PuV)Pt) = 1 9(Pt,v)v2{pt,v)p 2 d 2 p t dr] 
{Pt) J 9(Pt,v)Ptd 2 Ptdr] 

Again, the indexes F and Q indicate that all integrations in Eqs. f fl5|) . ffl6|) are 
performed over full phase space (F) or the phase space where particles are measured 
(O), respectively For completeness let us add that denotes the total number 
of produced particles (charged and neutral) and (p*) n is the average transverse 
momentum of the measured particles. 

One important lesson from the above result is that even if we measure particles 
in a limited fraction of the full phase space, e.g. a narrow pseudorapidity bin, the 
effect of transverse momentum conservation on (cos(0i + fa)) is not suppressed^]. 
This may be roughly understood in the following way: if each particle is generated 
in an independent manner except the constraint from overall transverse momentum 
conservation, then for each given particle its p t effectively has a chance oll/N (in 
the large iV limit) to be balanced by every other particle. That is equivalent to say 
that each pair has a back-to-back correlation of strength 1/N, which is preserved 
despite what fraction of particles is selected for measurement. Furthermore, while 
not changing the order of magnitude of the effect, the details of the p t and rj de- 
pendence of single particle distribution and v 2 may slightly affect the quantitative 
results. 

Next we examine two approximations to the main result, Eq. (|14p . 



5 Here and in the following we assume that azimuthal part of fl covers full 2tt i.e., <fi G [0, 2%). 
This is the only restriction we impose on f2. 

6 This statement is of course limited by the applicability of the central limit theorem, which 
however works for as little as 10 particles. 
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(i) If all the produced particles are measured, i.e. Q = F, and all have the same 
magnitude of the transverse momentum i.e. g(pt,v) ^ $(Pt ~ Po)h(r]), we 
recover the result (jSJ) discussed in the Introduction. 

(ii) If one allows for a finite acceptance but neglects the pt and r] dependence of 
t> 2 , i.e. v 2 ,f/q = v 2 ,f/q = v 2 ,f/q, Eq. (JUj) reduces to 

(cos(0 1 + 2 )) = -^Mi. (17) 

/ 1 2 

Since depends only weakly on the acceptance, Q, the corrections due to 

transverse momentum conservation already pointed out, more or less 

independent of the number of observed particles. 

Next let us calculate the contribution from transverse momentum conservation 
to the correlation function 5 = (cos(fa — 4>2)) which has also been measured by 
STAR. Following similar procedures, we obtain the following result: 

(cos(0i - fa)) = -- ' 2 ' '-, (18) 

We notice a few interesting features. First, the correlation scales like —1/N as 
expected. Second, the effect does not depend on elliptic flow in leading order, and 
thus, is much stronger than that in (cos(0i + fa)) which is of order Oiv^/N). Let 
us again examine the same two approximations to the above result. 

(i) If all the produced particles are measured, i.e. Q = F, and all have the same 
magnitude of transverse momentum i.e. g{pt,v) ^ ${Pt ~ Po)h>{v)> we obtain 



(cos(0i - fa)) 



'jv 



(ii) If one allows for a finite acceptance but neglects the p t and r] dependence of 
v 2 , i.e. ^f/u = V2,F/n = v 2 ,F/a, we obtain 

(cos^-^H-lMi. (19) 

Thus the main difference between the two correlators, 7 (Eq. IT4l) and 5 (Eq. 
T8|) is the presence of v 2 in the former. Therefore, effects of transverse momentum 
conservation will be much more visible in the correlator (cos(fa — 02)}- 



8 



3 Comparison with data 



In this Section we compare our results (1T4"|) and (JTSJ) with recently published STAR 
data [7] , where charged particles have been measured in the pseudorapidity interval 
— 1 < i] < 1 and the transverse momentum region p t > 0.15 GeV. For the results 
integrated over transverse momentum additional cut was imposed p t < 2 GeV. As 
seen from Eqs. (I14|I15|I16|[T8~]) to calculate the contribution of the transverse mo- 
mentum conservation to (cos(fa + fa)) and (cos(fa — fa)) we need full information 
about single particle distribution g(pt,v) an d elliptic flow v 2 (pt, if) in the full phase 
space. Unfortunately such complete information is not currently available. Thus 
we will make some reasonable assumptions that hopefully allow us to obtain an 
approximate insight into the discussed effect. 

First let us estimate the total number of produced particles N. From the PHO- 
BOS measurement (25] we know that the total number of charged particles N c h 
grows linearl}0 with the number of participants N part (or equivalently number of 
wounded nucleons [26J). At ^/s NN = 200 GeV we find that N ch « 14N part [23], thus 
the total number of particles can be reasonably approximated by 

N « (3/2)N ch = 21N part . (20) 

This result allows us to roughly estimate the contribution of transverse momen- 
tum conservation. We simply assume that {pf)ci = (pDf anc ^ ^2,n — ^2,f which leads 
to 

(cos(fa + fa)) ■ N part = -v 2 , n /21 w -0.004, (21) 

and 

(cos(fa - fa)} ■ N part = 1/21 w -0.05, (22) 

where we take W2,q = 0.08 to be slightly larger than the elliptic flow parameter 
V2,n = 0.06 in order to account for the momentum dependence of v%. Below we 
will show that these simple assumptions are well reproduced in a more detailed 
calculation. 

While this effect gives a contribution with the same sign and order-of-magnitude 
for the same-sign pair correlation data, it is a factor 3 — 5 (very peripheral - mid- 
central) less in magnitude for (cos(fa + fa)) and a factor 1.5 — 4 (mid-central - very 
peripheral) larger than the STAR data. It also gives the same contribution to the 
opposite-sign pair correlation for which the data are positive. 

To perform more precise calculations we further assume that the single particle 
7 In contrast to the number of charged particles at midrapidity that grows slightly faster than 

N part • 
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distribution g(pt, rf) can be expressed in the following way 



g(p t ,v) = 7^ ex P ( -TFT ) Kv), (23) 



T 2 \ T 

where for simplicity we take a thermal distribution and assume factorization of the 
momentum and pseudorapidity dependence. To introduce a slight pseudorapidity 
dependence of T [27j we will distinguish between T F and T mid depending upon 
whether we integrate over full pseudorapidity range or the midrapidity region. In 
the following we take T mid = 0.225 GeV and T mid /T F = 1.1 [27]. In Eq. ([23} 
h(i]) is the normalized pseudorapidity single particle distribution, which can be well 
represented by a double Gaussian form 



h{rj) 



1 



2V27ra 



2 



exp — + exp 



(24) 



2a 2 J r V 2a 2 

with r] = 2 and o = 1.9 for y/s^ = 200 GeV@ Finally we have to specify the mo- 
mentum and pseudorapidity dependence of elliptic flow. For simplicity we represent 
this dependence with a factorized liner ansatz for both the transverse momentum 
and the pseudorapidity, which represents the presently available date reasonably 
well up to p t ~ 2 GeV [291 M M\ B 

7 — \tj\ 

MPt,v) = Cp t — - — , (25) 

and t>2 = for > 7. Here C = 0.14 so that for the midrapidity region \r)\ < 1 the 
calculated elliptic flow equals 0.06 [291 EH E] ■ 

With these parametrisations we find that: (pt)1i ~ 0.257, (p 2 ) F ~ 0.252, v^n ~ 
0.0886 and v 2i f ~ 0.0773. Substituting these number into Eqs. f[T4"|) . f lTH]) we obtain 

(008(0! + 2 )) • iV part w -0.005, (26) 

and 

(cos(0! - 2 )> • N part w -0.05, (27) 

which are very close to the numbers estimated at the beginning of this Section. 
While it is difficult to estimate the precise uncertainty of our calculation, however, 
we expect our results to be correct within a few tens of percent. 



4 Differential distributions 

The STAR collaboration also presented data for (cos(0i + 2 )) as a function of 
p + = (pi tt + p2,t)/2 and p_ = \pi tt — P2,t\- Both distributions are very informative: 

8 We performed the fit to the PHOBOS 35 - 45% centrality AuAu data [28] . 

9 We have checked the influence of constant for p t > 2 GeV and found it negligible. 
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the azimuthal correlation (cos(0i + 2 )) increases roughly linearly with p + while it 
depends only weakly on p_. In [JJ] we showed that such behavior is not inconsistent 
with the CME and can be understood if we assume that correlated pairs have slightly 
larger momenta than uncorrelated ones. 

As we will show below, a qualitatively similar behaviour can be obtained from 
transverse momentum conservation. With the two-particle distribution ( ITT]) , the p + 
differential distribution reads 

/ fc{Pl,P2)S (2p+ - Pl,f - P2,t) COS(01 + 2 )^ 3 Pl^ 3 P2 
(008(0! + 2 )) p+ = ^ ^ ( 28 ) 

n 

and analogously for p_ = — p 2 ,t|- 

Let us first discus the simplified case where V2{pt,v) is replaced by its average 
value i>2- Taking Eqs. ffTTj) and f )23|) into account we obtain^] 

2d 2 

(cos(0 1 + 2 )) p+ = -t; 2I ^, (29) 

and 

(C0B(& + 2 )} = 3T2 + 3T P-+P- (30) 

where for simplicity we set Tp = T m id = T. As can be seen, the correlation shows 
a strong quadratic growth with increasing p + , in qualitative agreement with data. 
The dependence on p_, on the other hand, is essentially constant for p_ <T before 
it exhibits a linear increase. Performing analogous calculations for S we obtain 

(COS(01 - 02)) p+/p _ = ^ (COS(01 + 2 )) p+/p _ , (31) 

thus the p + and p~ dependence is identical however, the signal is significantly 
stronger. 

We also performed full calculations with g(pt,i]), h(r)) and f 2 (pt,^) discussed in 
the previous Section and given by Eqs. ( |23l) . ( 1241) and ( 1251) . respectively. However, in 
this calculation we correct v 2 (pt, v) an d assume a constant value for v 2 for transverse 
momenta p t > 2 GeV. The results for N part = 100 are presented in Fig. CD and Fig. 
[2] for (cos((/>i + 2 )) and (cos(0i — 2 )), respectively. 

Again, similarly to what has been observed in the STAR data, the correlation 
(cos(0i + 2 )) grows rapidly with increasing p + while appears much flatter with 
increasing p_. Comparing our results with the STAR data we see that the transverse 



10 In this section for simplicity we integrate from pt — despite the finite cut pt > 0.15 GeV in 
the STAR experiments. 
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0.0 0.5 1.0 1.5 2.0 

(Pi,f + 02,f)/2; IPi.f-P2.fi (GeV) 

Figure 1: (Color online) The two-particle azimuthal correlation (cos(0i + fa)) vs 
P+ = (Pi,t +P2,t)/2 (blue line) and p_ = \p ljt - p 2 ,t\ (red line) for N part = 100. 

momentum conservation gives comparable signal for p + and p- larger than 1 GeV 
and underestimate the data for lower values of p + and p-. 

Finally let us discus the pseudorapidity dependence of the correlation (cos(0i + fa)). 
The dependence of {cos(fa + fa)) on I771 — 772 1 has been measured by STAR and 
found to be dominated by \rji — 772 1 < 2, which is consistent with the CME expec- 
tation. It is quite clear that such dependence cannot be obtained in the present 
calculation, as no pseudorapidity dependence appears in the nontrivial part of the 
two-particle correlation function shown in Eq. flTTJ). Consequently, transverse mo- 
mentum conservation predicts the correlator (cos(fa + fa)) to be essentially flat as a 
function of | rji — 772 1 in the midrapidity region except for a very mild dependence due 
to the slight dependence of f(p) and v-i on 77. However, here we have assumed that 
the transverse momentum is balanced over the entire rapidity interval. In an actual 
heavy ion reaction it is not unreasonable to expect that the transverse momentum is 
balanced over a shorter rapidity interval. If this were the case, we would predict not 
only a stronger rapidity dependence of the signal but also a considerably stronger 
signal at midrapidity. Therefore, it would be worthwhile to construct and measure 
an equivalent of the charge balance function [32] for the transverse momentum. This 
problem is currently under our consideration. 

5 Conclusions and comments 

We have quantitatively investigated the contribution of transverse momentum con- 
servation to the azimuthal correlation observables 7 = (cos(0i + fa)) and 5 = 
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0.0 0.5 1.0 1.5 2.0 

(Pi,f + P2,f)/2; IPi,f-p 2 ,fl (GeV) 

Figure 2: (Color online) The two-particle azimuthal correlation (cos(0i — 02)) vs 
P+ = (pi,t+P2,t)/2 (blue line) and p_ = -p 2 ,t| (red line) for i\r port = 100. 

(cos(0i — 02)) measured by the STAR collaboration as motivated by the possible 
strong local parity violation and Chiral Magnetic Effect. Our conclusions can be 
summarized as follows. 

(i) The contribution due to transverse momentum conservation is comparable 
in magnitude to the prediction of the Chiral Magnetic Effect as well as the 
data. In the STAR acceptance we find this contribution to be approximately 
equal to 7 ~ —l.lv^/N, where t> 2 is the elliptic flow coefficient at midrapidity 
and N is the total number of produced particles (neutral and charged). This 
result suggests rather week energy dependence at RHIC since v 2 and N scale 
similarly with energy. 

(ii) Taking v 2 = 0.06 and N = 21N part where N part is the number of participants 
we obtained 7 • N part « —0.005, which is a factor 3 — 5 (very peripheral - mid- 
central) smaller than the experimental data. Thus we may conclude that the 
transverse momentum conservation alone cannot explain the data. Also there 
is no charge-dependence as opposed to experimental data. It is however a 
significant source of background that eventually must be quantified and taken 
into account if one really wants to extract the possible CME from the present 
(and future) data. 

(iii) We have demonstrated that finite acceptance issues, i.e. the facts that particles 
are measured in a relatively narrow pseudorapidity bin and neutral particles 
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are not detected, do not suppress the effect of transverse momentum conser- 
vation on 7. 

(iv) We studied the dependence of 7 vs p + = (pi ;t +p 2 ,t)/2 and p_ = \p\ tt — p 2 ,t|- 
We found that 7 increases with increasing p + , 7 oc (^y) a+ with a + = 2 — 3. 
The dependence onp_ is much weaker 7 oc (^y) a ~ with a_ ~ 1. This behavior 
is qualitatively similar to what is observed in the data. We also investigated 
the dependence of 7 on [771 — 77 2 | and found no pseudorapidity dependence in 
contrast to what is observed in the data. 

(v) Finally we calculated 5 = (cos(</>i — 2 )) and found that in the STAR accep- 
tance 5 ~ —1/N. We found this contribution to be a factor 1.5—4 (mid-central 
- very peripheral) larger than the experimental data indicating again that the 
transverse momentum conservation effect is a significant source of background. 

(vi) The present calculation is based on the minimal assumption that the transverse 
momentum is balanced over all particles (in the full phase space). Thus it is 
very likely that the calculated contribution to 7 and 5 from the transverse 
momentum conservation represents rather the lower limit. Should, as it is 
not reasonable to assume, the transverse momentum be balanced over a finite 
rapidity interval, we predict not only a stronger effect at midrapidity but also a 
rapidity dependence of the correlation function 7 and 5. Thus a measurement 
of something like a transverse momentum balance function would be highly 
desirable. This problem is currently under our consideration. 

We end this paper with a few further comments. 

(a) While transverse momentum conservation alone is not sufficient to explaining 
the data, one may combine it with other effects such as the Chiral Magnetic 
Effect or Local Charge Conservation [10] to get closer to the data. In Table. [1] 
we summarize the estimated contributions to the azimuthal correlations from 
these effects together with the STAR data. All numbers quoted are for AuAu 
200 GeV collisions at about 50 — 60% centrality (corresponding to N part « 50 

As a precaution, all numbers (except the STAR data) bear considerable un- 
certainty. The numbers from Chiral Magnetic Effect (CME) for observables 
(cos(0i + 4>2)) ++ + _ are extracted from [4J. The numbers for (cos(0i — </> 2 )} ++ + _ 
where obtained using the relation (cos(0i — <f> 2 )} ++ = — (cos(0i + 2 )) ++ and 
analogously for (+, — ), which hold in case of a pure CME. The numbers from 
Local Charge Conservation (LCC) are inferred from [TU]: the authors showed 
the difference (cos(0i + 2 )) + _-§ [(cos(0i + 2 )) ++ + (cos(0i + 2 )) ] • From 
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x 10 3 


(cosOi + 2 )) ++ 


(COS(01 + 4> 2 )) + _ 


(cos(0i - 4> 2 )) ++ 


(COS(0! - <p 2 )) + _ 


CME 


-(0.1-1) 


+ (0.01-0.1) 


+ (0.1-1) 


-(0.01-0.1) 


LCC 


~ 


+ (0.1-1) 


~ 


+ (1-10) 


TMC 


0.1 


0.1 


— 1 


~ -1 


DATA 


-0.45 


+0.06 


-0.38 


+1.97 



Table 1: Estimated contributions to azimuthal correlations from various effects and 
comparison with data. CME, LCC, and TMC stand for Chiral Magnetic Effect, 
Local Charge Conservation, and Transverse Momentum Conservation, respectively, 
while the DATA is from STAR measurement for AuAu 200 GeV collisions at about 
50 - 60% centrality. 



LCC one should expect (cos(0i + <t> 2 )) ++ ~ so the (cos(0i + <t> 2 )) + _ 

could be inferred. Furthermore the value for (cos(0i — <p 2 )) + _ is estimated by 
(cos(0i — <t> 2 )) + _ ~ (cos(0i + 4> 2 )) + _ /v 2 . Finally as the authors pointed out, 
their results represent an upper limit of the magnitude for LCC as they enforce 
exactly local charge neutrality. Our results for Transverse Momentum Conser- 
vation (TMC) are as given in the previous sections (with the expected uncer- 
tainty within a few tens of percent). The STAR data is from [7J. Even though 
these are rough estimates, one still can make a few observations: first, no single 
effect shows a pattern for all observables that are in accord with data; second, 
different correlators appear to be dominated by different effects, in particular 
- both CME and TMC provide important contributions to (cos(0i + </> 2 )) ++ , 
while TMC seems necessary for explaining the (cos(0i — <fr 2 )) ++ , and LCC for 
the observed value of (cos(0i — <p 2 )) + _. The situation for (cos(0i + 2 )) + _ on 
the other is more complicated and none of the effects discussed here seems to 
dominate. Clearly, additional measurements will be required to disentangle 
this situation. 

(b) As far as comparisons of the data with models is concerned, such as the ones 
presented by the STAR collaboration [7J, one has to ensure that the model 
satisfies at least two criteria. First, the model has to conserve transverse 
momentum not only on average but event-by-event. Second, the model has 
to reproduce the measured magnitude of the elliptic flow, v 2 , as essentially 
all "trivial" contributions to (cos(0i + </> 2 )} +++ _ scale with v 2 . For example, 
UrQMD is known to underestimate the measured v 2 [31]. This may be partly 
the reason that it also underestimates the measured data for (cos(0i + <p 2 )) ++ 
as reported in Ref. [7J. 
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